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ABSTRACT 

A properly calibrated longitudinal magnetograph is an instrument that measures circular polarization and gives an estimation of the 
magnetic flux density in each observed resolution element. This usually constitutes a lower bound of the field strength in the resolution 
element, given that it can be made arbitrarily large as long as it occupies a proportionally smaller area of the resolution element and/or 
becomes more transversal to the observer and still produce the same magnetic signal. Yet, we know that arbitrarily stronger fields are 
less likely -hG fields are more probable than kG fields, with fields above several kG virtually absent- and we may even have partial 
information about its angular distribution. Based on a set of sensible considerations, we derive simple formulae based on a Bayesian 
analysis to give an improved estimation of the magnetic field strength for magnetographs. 
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1. Introduction 


The activity of the Sun (and of a large percentage of other stars) 
is driven by the presence of magnetic fields. It was realized long 
ago that it is important to carry out systematic observations of 
the surface magnetism of the Sun if we ever want to understand 
its activity cycle and that of other stars. This possibility become 
a reality after the discovery of the Zeeman effect ( |Zeeman|1897] ) 
that a magnetic field produces on some spectral lines (leading 
to a splitting in frequency of the <x and n components), and its 
discovery on the Sun by [Hale| (| 1908 i 


A breakthrough occured when 


Babcock ( 1953| > developed 


the longitudinal magnetograph, a device that measured the line- 
of-sight (LOS) component of the magnetic field over an area of 
the surface of the Sun. In order to measure the magnetic field, 
such apparatus performs a measurement of the intensity 1(A) and 
circular polarization V(A) of a magnetically sensitive spectral 
line in a certain narrow filter with transmission p(A), and returns 
the following quantity (e.g., |Landi DeglTnnocenti & Landolfi| 
2004) : 


S v - 


f dAV(A)p(A) 
f d AI(A)p(A) ' 


(1) 


The importance of this quantity resides on the fact that it can be 
easily related to the magnetic flux density (<\>) across the reso¬ 
lution element of the instrument (pixel Q To this end, a number 
of simplifying conditions have to be imposed, the most impor¬ 
tant of which is the assumption that the field is in the weak-field 
approximation (e.g., |Landi DeglTnnocenti & Landolfi|2004| ). If 
this holds, the measured signal is simply given by 

Sy = Cfi>, (2) 

where C is a calibration constant that needs to be obtained from 
a proper modeling of the line formation mechanism and depends 


1 Although the units Mx cm -2 and G are equivalent, we follow Keller 
|et al.| ( |1994| > and use the former for the magnetic flux density to make 
its observational character explicit. 


on the exact details of the filter p(A) and the selected spectral 
line. Once the instrument is calibrated, fi) = S y/C gives the esti¬ 
mated magnetic flux density. When the calibration is assumed to 
be perfectly known and the noise in S y is Gaussian, <I> = S y/C 
constitutes the unbiased maximum likelihood estimation of the 
magnetic flux density. 

From its conception, many longitudinal magnetographs have 
been on routine operation around the globe and on space. One 
of the most used is the Michelson Doppler Imaging (MDI) on¬ 
board the Solar and Heliospheric Observatory (SoHO), which 
provided synoptic maps of magnetic flux density. Obviously, the 
fundamental drawback of any longitudinal magnetograph is that 
it only gives us partial information about the magnetic field in 
the observed region (Lites et al.] |1999[ ), giving a lower limit to 
the strength. We propose in this note, using simple Bayesian in¬ 
ference, a way to probabilistically correct the measurements of 
longitudinal magnetographs to provide a better estimation of the 
magnetic field strength, instead of just a lower limit. The correc¬ 
tion relies on some sensible statistical a-priori assumptions about 
the properties of the field in the observed resolution element that 
can be modified at will. This facilitates the adaptation of our sug¬ 
gested correction to any desired a priori assumption imposed by 
the researcher. 


2. Bayesian analysis 

2.1. Generative models 


The first step in any Bayesian analysis is to write down an ex¬ 
pression that describes how the observations are generated, the 
so-called generative model. If the field in the observed resolu¬ 
tion element is resolved at the spatial resolution of the instru¬ 
ment (i.e., the magnetic field is assumed to be homogeneous and 
filling the whole resolution element), then the magnetic flux den¬ 
sity can be directly related to the magnetic field strength B in the 
weak-field limit (e.g., Landi Degl’Innocenti & Landolfi|2004| >. 
Without loss of generality, we assume that C = 1, so that our 
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measurement is O = Sy. This avoids carrying the instrument- 
specific constant C in all calculations. Therefore, the generative 
model for the resolved case is 

eh,- - Bp + e, (3) 

where p = cos 6 is the cosine of the angle between the LOS and 
the magnetic field vector, while e stands for the measurement 
noise, which we assume to have Gaussian statistics with variance 
cr 2 . Only measuring <t> r , a degeneracy between B and p appears 
that cannot be resolved unless more information is added. We 
show in the following that the inclusion of priors helps partially 
resolving this degeneracy. 

When the field in the resolution element of our instrument 
is not resolved, it is useful to introduce the concept of a filling 
factor, / e [0,1], This quantity accounts for the fraction of the 
resolution element that is filled with an homogeneous magnetic 
field. The remaining 1 - / fraction is assumed to be field-free. 
Under the assumption that the presence of a magnetic field does 
not modify the thermodynamic conditions, the generative model 
for the non-resolved case is 

<h nl . = fBp + e. (4) 

It is possible to generalize these formulae for the case in which 
the thermodynamic conditions are different in the magnetic and 
the non-magnetic regions ( |Landi Degl’Tnnocenti & Landolfij 
|2Q04) >, but we prefer to stick to these simple cases. Our analy¬ 
sis is fully generalizable to these cases, though. 


2.2. Posterior distribution 


Once the magnetic flux density fiJ obs is measured, all the infor¬ 
mation about the model parameters of interest is encoded in the 
posterior distribution. The posterior distribution for the genera¬ 
tive model of Eq. 0 is, applying the Bayes rule, given by 


P(B,p,f |®obs) 


p(<t> obs \B,p,f)p(B,p,f) 

P(® obs) 


( 5 ) 


where p(B,p,f |<t> 0 b s ) is the posterior for the model parameters, 
p(®obs\B,p,f) is the likelihood that encodes the ability of the 
model to reproduce the observations, p(B, p. f) is the prior dis¬ 
tribution that summarizes all the a-priori information we have 
about the parameters. Finally, p(0 0 bs) is the evidence, a normal¬ 
ization constant that is unimportant here because it does not de¬ 
pend on the model parameters and that we drop in the following. 

According to the statistics of the noise proposed in the gen¬ 
erative model, the likelihood is given by a Gaussian distribution 
with variance cr 2 : 


Z = p(<t>obs\B,p,f) = —=— exp 

y2n<j n 


(O obs - Bfp ) 2 

2cr2 


( 6 ) 


For the sake of simplicity, we assume that B, p and / are a-priori 
statistically independent, so that the prior factorizes: 


P(B,p,f) = p(B)p(p)p(f). (7) 

Note that this step is just a simplification and more complicated 
a-priori information can be easily introduced. In particular, it is 
possible to easily introduce the fact that stronger fields tend to 
be more vertical by writing a combined prior p(B,p ) that takes 
this property into account. 

Our state of knowledge about the magnetic field strength af¬ 
ter carrying out the observation with the magnetograph is com¬ 
puted by integrating / and p from the posterior distribution of 
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Eq. (|5j. Assuming complete ignorance on / by using a flat prior 
on the interval [ 0 , 1 ], the ensuing marginal posterior is: 

p(B\<P ohs ) oc p(B) J dfp(f) J dpp(p)p(<$ ohs \B, p, f) 

= P(B ) f d ^)^-[erf(^)-erf(^^ 
J-1 Bp \ V2 cr,,) \ V2cr„ 


( 8 ) 


Once the posterior distribution is computed, it is easy to estimate 
upper limits to the magnetic field strength at a certain credibility 
limit. We later calculate Bgs and B 95 , as the 68 % and 95% upper 
limits. 

It is key to understand that the posterior distribution for the 
magnetic field strength will change depending on the election of 
the specific prior. Although this is obviously always the case, it 
is of special relevance in our problem because of the presence of 
strong degeneracies. For this reason, it is advisable to consider 
several priors and see how the results change with them. The first 
prior we consider is a log-normal distribution, as pointed out by 
Vogler (2003|> and Vogler et al. (20051) using magnetoconvection 


simulations: 


P(B) 


Bcr b V2^r 


exp 


(In B- In B 0 ) 


2cr 2 b 


(9) 


where Bq and <x/, are free hyperparameters, that we choose to be 
Bo = 38 G and cr/, = 1.2, following the results of (Dominguez 
Cerdena et a l.|2006j >. Note that the average of a log-normal ran¬ 
dom variable equals exp(ln B {l + cr 2 /!), which equals 78 G for 
the values of the parameters. Although this value is slightly be¬ 
low those estimated by Trujillo Bueno et al.| ( |2004] l, |Shchukina| 
|& Tru jillo Bue no| ( |201 1 1 and Rempel ( 2014[ >, it decreases slower 
towards large magnetic field strengths than the exponential dis¬ 
tribution present in the simulations or assumed by Trujillo Bueno 
|et al.|(|2004|. Our second prior is an exponential distribution 




( 10 ) 


with Bo = 130 G ( Trujillo Bueno et al.|2004| Rempel|2014| l or 
Bq = 60 G (Rempe l|2014| i. We advance that, although the three 
considered priors are very different, the results are relatively ro¬ 
bust to this election. 


3. Results 

3. 1. Ignorance about the topology of the field 

When we are completely ignorant about the topology of the field, 
a reasonable assumption is to consider it, a-priori, isotropic. To 
this end, we consider a flat distribution in p, which yields a vec¬ 
tor field uniformly distributed in solid angle. Figure [T] displays 
the marginal posteriors for a fixed standard deviation of the noise 
of cr„ = 5 Mx cm 2 , typical of modem longitudinal magne¬ 
tographs 2 ] We show the marginal posteriors for different values 
of the observed magnetic flux density. The left panels display 
the marginal posterior, while the right panels show the cumula¬ 
tive distribution. The computation is done using a proper numer¬ 
ical quadrature in p. Each row corresponds to different priors, as 
discussed above. We note that the effect of the specific prior is 

2 The code to reproduce the figures can be found in 
https://github.com/aasensio/magnetographCorrection. 
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Fig. 1 . Marginal posterior distributions from Eq. {8} (left panel) and their ensuing cumulative distributions (right panel) for a standard deviation 
of the noise of cr„ = 5 Mx cnr 2 and different values of the observed magnetic flux density. The actual values of ® obs are indicated in the legend. 
Flat priors for / and yu are used. Each row shows the calculations for the different priors: (top) a log-normal prior with So = 38 G and oy, = 1.2, 
(middle) an exponential prior with B 0 = 130 G, and (bottom) an exponential prior with B 0 = 60 G. 


somehow marginal, except in the tails of the posteriors. It is ob¬ 
vious that the exponential prior with Bq = 60 G disfavors strong 
fields and the tails of the posterior become less important. 

The curves of Fig. [T] are conditioned on the actual measured 
value of ® 0 b s . Consequently, consecutive measurements of ex¬ 
actly the same resolution element would lead to slightly differ¬ 
ent curves. In order to see the influence of the noise, we show in 
Fig.Elthe marginal posteriors for a fixed value of <F ohs = 40 Mx 
cm ~ and for different values of cr„ for the three considered pri¬ 


ors. The figure demonstrates that, when the noise level is small, 
values of B below <I* ohs have very small probability. When the 
standard deviation of the noise increases, values of B < ( hobs 
gain some probability. This result clearly exposes that, when the 
standard deviation of the noise is of the order of the measured 
magnetic flux density, <E> 0 b s is not anymore a lower bound of the 
magnetic field strength. The differences among the priors are not 
very large except in the tails. This shows the robustness of the 
inference to the specific chosen prior. 
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Fig. 2. Same as Fig.[l]for a fixed observed magnetic flux of O 0 bs = 40 Mx cm 2 and different values of the standard deviation of the noise <x„. 


The marginal posteriors indicate that there is some available 
information about the magnetic field strength on a single mea¬ 
sure of the magnetic flux density. In essence, small values of 
B are preferred with respect to large values. The peak of the 
marginal posterior p(B |<f> 0 t, s ) (the so-called marginal maximum 
a-posteriori [MMAP] value) for the magnetic field is shifted to¬ 
wards increasingly larger values when the measured <I* ohs in¬ 
creases. We extracted the MMAP values of the magnetic field 
strength from these curves and we show them in the upper pan¬ 
els of Figure [3] with respect to the measured magnetic flux den¬ 
sity for different values of the noise and for the three considered 
priors. For the cases with standard deviations of the noise be- 
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low 10 Mx cm~ 2 , it follows a relatively linear trend. For the re¬ 
maining ones, large variations are a consequence of the specific 
noise realization in the displayed posterior (it will change from 
realization to realization). We have displayed a power-law fit to 
the curves. Our proposal is to use this value as an estimation of 
the most probable value of the magnetic field strength when the 
effect of the inclination and filling factor are considered. This 
MMAP value is very robust to the specific hyperparameters of 
the prior p(B). 

A problem that arises with the MMAP value is that it is bi¬ 
ased towards small values of the field strength, given that the 
posterior distribution is heavily skewed. More information can 
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Fig. 3. Marginal maximum a-posteriori (upper panel), percentile 68 (middle panel) and 95 (lower panel) obtained from the marginal posteriors 
computed from Eq. ©■ Each curve corresponds to a different value of the standard deviation of the noise. The dashed lines are the fits to the 
curves, a straight line for Bmmap and a function of the type o-O^ for B 68 and Bg 5 . 




Fig. 4. Same as Fig. Qfor ^ = 5 Mx cm 2 and ® obs = 40 Mx cm 2 and different values of the anisotropic index of the prior for y j. 


be extracted, though, from the cumulative distribution functions, 
shown in the right panels of Figs. [1] and [2] because they allow us 
to put strict upper limits to the magnetic field. From these curves, 
we have extracted the percentiles 68 and 95 and they are shown 
in the second and third panel of Fig. [3] Again, when the noise is 
not too large, they approximately fulfill the power-laws shown in 
each panel. These percentiles are also quite robust to the specific 
election of the prior and its hyperparameters, although slighly 
more sensitive than the MMAP value. 

3.2. Anisotropic field distributions 

It has been argued in the literature, both observationally (Sten- 
|flo|2010| > and from simulations (Schiissler & Vdgl er|2008} , that 


the distribution of fields might be anisotropic. A simple way of 
imposing this a-priori is by modifying the flat prior for /./ to 

p(M)*\p\ a , (ID 

where the parameter a controls the anisotropy (a — 0 is an 
isotropic distribution, a —» 00 corresponds to vectors aligned 
with the LOS, while a < 0 force fields to be predominantly or¬ 
thogonal to the LOS). Figure [4] displays the marginal posterior 
for different values of a for the log-normal prior. For increasing 
values of a, the tails of the marginal posterior become less and 
less important, converging in the limit a —> 00 to a Gaussian dis¬ 
tribution with mean <t> 0 b s and variance <x 2 , i.e., what one would 
obtain by carrying out the inference assuming that |/r| = 1. For 
not very anisotropic distributions, the MMAP value is very ro- 
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Fig. 5. Posterior distribution for the transversal component of the field for different values of ®„b s and cr„. The value of the upper limit transverse 
field at 95% credibility is shown in each panel. The bin size is optimally chosen for each case. 


bust to the specific anisotropy of the field distribution, so that the 
corrections derived above can be reliably applied. 

3.3. Transversal component 

The Bayesian analysis we have carried out also allows us to 
put some constraints on the transversal component of the mag¬ 
netic field, B ± . This might look surprising but it is obvious that 
putting limits to B automatically puts limits on B. Give that 
the transversal component is given as a change of variables by 
IT = Bfyfl - p 2 , it is not possible to do much advance ana- 
litically. For this reason, we consider a numerical solution to 
the problem by sampling from the full posterior of Eq. 0 US- 
ing a Markov Chain Monte Carlo method and compute the 
marginal posterior from the samples. The results are shown in 
Fig. § where we consider different values of the observed mag¬ 
netic flux density and standard deviation of the noise. Each panel 
shows also the estimated value of the transversal component of 
the field at 95% credibility. 

It is obvious that the correct way of having an estimation of 
the transversal component of the field is by observing the lin¬ 
ear polarization components of the Stokes parameter. However, 
the results of Fig. [5] provide an upper limit based on sensible 
assumptions and on the observation of the magnetic flux den¬ 
sity alone. When linear polarization is available, it is necessary 
to augment the generative model to include this information, as 
done by Asensio Ramos (|2011 ). 


3 We use the emcee package developed by Foreman-Mackey et al. 
< 20l3) . 


3.4. Application to real data 

We apply the previous formalism to the data obtained with IMaX 
(Martinez Pill et et aT7| |201 1[ > onboard Sunrise ( |Solanki et al.| 
12010) on the first flight on June 2009. Although IMaX is a vec¬ 
torial magnetograph, that measures the four Stokes parameters, 
we only make use of the inferred magnetic flux density ® Q b s in 
our estimation of the magnetic field strength. One of the partic¬ 
ularities of the IMaX mission is that both the original maps and 
those reconstructed with the phase-diversity post-facto recon¬ 
struction algorithm were provided (Paxman et al.|[T992| [Vargas 
Dommguez|2009). The left column of Fig.[6]displays the maps of 


4>obs in the two cases, demonstrating that the magnetic structures 
in the reconstructed maps are more compact than in the original 
ones. To have a rough estimation of the magnetic field strength, 
we compute it as (Bj + /J 2 ) l/2 , where B\\ and B_ are obtained us¬ 
ing the simple calibration curves of Martinez Pillet et al. ( 2011) . 
This is an imprecise but probably conservative approximation 
for a filter-based instrument. This estimation is displayed in the 
middle column. As noted by Martinez Gonzalez et al. ( |2012) , the 
maximum-likelihood estimation of B ± is biased in the presence 
of noise (because B ± is obtained from ( Q 2 + f/ 2 ) 1 ^ 2 , a positive 
definite quantity even in the absence of signal), which results in 
the appearance of a pattern that mimics inverse granulation. Fi¬ 
nally, the right columns of Fig. [6] display the upper limit at 95% 
credibility using our estimations for an exponential prior with 
B{)- 130 G and a very small noise variance (see the lowest cen¬ 
tral panel of Fig. [3}. The upper limit field is larger for the decon¬ 
volved data, a consequence of correcting for the spatial spread 
introduced by the instrument. 
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Fig. 6. Original (upper panels) and deconvolved (lower panels) inferred parameters for the IMaX data of the first flight. The first column shows 
the magnetic flux density. The middle column displays the inferred magnetic field strength, obtained as (Bjj + B 2 ) 1/2 . The last column displays our 
inferred upper limit to the field strength at 95% credibility. 
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Fig. 7. Comparison of the histogram (upper panel) and cumulative his¬ 
togram (lower panel) of the magnetic field strength computed using a 
standard inversion code and our Bayesian approach at 95% credibility. 


As a further check of our approach, we have compared the 
upper limit at 95% credibility for the magnetic field strength with 
the results of a standard inversion with the SIR code (Stokes In¬ 
version based on Response function; |Ruiz Cobo & del Toro Ini-| 
|esta| 1992]) of the data (kindly provided by L. Bellot Rubio) using 
a simple model comparable to a Milne-Eddington atmosphere. 
Figure [ 7 ] displays the results in the form of histograms (upper 
panel) and cumulative histograms (lower panel). The results are 
very similar, even reproducing the position of the peak. 


4. Conclusions 

We have shown that the application of Bayesian inference to the 
data obtained with longitudinal magnetographs can give some 
information about the magnetic field strength. This information 
is a product of the presence of priors on the magnetic filling fac¬ 
tor and the inclination of the field. Our calculations have been 
done with very uninformative priors on the angular distribution 
of the magnetic field and prior distributions for the magnetic field 
strength extracted from analysis of simulations and observations 
of the quiet Sun. We demonstrate that our results are relatively 
robust to the specific election of the prior. 

The calculations can be repeated introducing more informa¬ 
tive priors for /, // and B in Eq. (|HJ. Irrespectively of the chosen 
prior, Eq. ([8]» is easy to compute using any suitable numerical 
quadrature for two-dimensional integrals. 
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